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Motivated by the anomalous positron flux recently reported by the PAMELA collaboration, wc 
study the cosmic-ray positron produced by the pair annihilation and the decay of superparticle dark 
matter. We calculate the cosmic-ray positron flux and discuss implications of the PAMELA data. 
We show that the positron excess observed by the PAMELA can be explained with some of the 
superparticle dark matter. 



Dark matter, which accounts for about 23 % of the 
mass density of the present universe [l[ , has been a mys- 
tery in the fields of particle physics, astrophysics, and cos- 
mology. Even though there exist various well-motivated 
candidates for dark matter, like the lightest superparti- 
cle (LSP), axion, and the lightest Kaluza-Klein particle 
in the universal extra-dimension scenario, the nature of 
dark matter is completely unknown and many possibili- 
ties have been discussed to study dark-matter properties 
with direct and indirect detections of dark-matter signals 
as well as with colliders. Among them, measurements of 
the anti-particle fluxes in the cosmic ray give important 
information about dark matter; annihilation or decay of 
some classes of dark matter particle produces energetic 
anti-particle, which may be observed in the cosmic ray. 

Very recently, the PAMELA collaboration has an- 
nounced the first result of the measurement of the 
positron flux for the energy range 1.5—100 GeV 0], which 
indicates an anomalous excess of the positron flux over 
the expected background. In particular, the PAMELA 
results show a significant increase of the positron frac- 
tion at the energy range 20 GeV < E < 100 GeV, even 
though the positron fraction is expected to decrease as 
energy increases [H, BJ ■ In addition, in the past measure- 
ment of the positron flux by HEAT 0] , similar excess was 
also pointed out. To understand the anomalous behavior 
of the positron flux, an unaccounted mechanism is neces- 
sary from particle-physics and/or astrophysical point of 
view. 

One important possibility of realizing enhanced 
cosmic-ray positron flux is the annihilation or the de- 
cay of dark matter particle @, 0|- (Recently, another 
possibility is also pointed out that the anomalous flux 
may be due to the positron emission from Pulsars 
In this study, we pursue a possibility that the anomalous 
cosmic-ray positron is from the pair annihilation or the 
decay of the dark matter particle. There are various can- 
didates for dark matter. Among those, in this study, we 
consider one of the most well-motivated candidates for 
dark matter, i.e., the LSP in models with low-energy su- 
persymmetry (SUSY). Conventionally, in the LSP dark 
matter scenario, it is assumed that the i?-parity is con- 
served and that the LSP is Bino-like lightest neutralino or 
the gravitino. Then, the cosmic-ray positron flux is neg- 
ligibly small. This is because, in such a case, the positron 
is produced by the pair-annihilation of the LSP, whose 



cross section is significantly suppressed in the cases of the 
Bino-like and gravitino dark matter. However, the above 
assumptions are not necessary in realizing the LSP dark 
matter scenario; in several well-motivated scenarios, one 
of the above assumptions are relaxed and an enhanced 
positron flux may be realized. 

In this letter, we consider the implication of the 
positron-flux observation of the PAMELA experiment to 
the SUSY model in which the LSP becomes dark mat- 
ter. We calculate the production rate of the cosmic- 
ray positrons in several scenarios. Then, we show that, 
even if we adopt the conventional estimation of the back- 
ground e ± fluxes, the positron excess observed by the 
PAMELA can be explained in some of the cases. 

Let us first summarize our procedure to calculate the 
positron flux <I> e + (as well as the electron flux <E> e - ). (For 
detail, see We solve the diffusion equation to take 
account of the effects of the propagation of the positron. 
The energy spectrum of the positron from dark matter 
f e + (E, r) evolves as [1] 

df e + d 



K{E)V z f e 



[b(E)f e+ ] + Q. 



(1) 



dt v ' Je dE 

The function K is expressed as K = K Eq cV [10l |. where 
Eq c y is the energy in units of GeV. In our numerical 
calculation, we use the following three sets of the model 
parameters, called MED, Ml, and M2 models, which are 
defined as (6, if [kpc 2 /Myr] , L[kpc]) = (0.70,0.0112,4) 
(MED), (0.46,0.0765,15) (Ml), and (0.55,0.00595,1) 
(M2), with R = 20 kpc for all models. Here, L and R are 
the half-height and the radius of the diffusion zone, re- 
spectively. The MED model is the best-fit to the boron- 
to-carbon ratio analysis, while the maximal and mini- 
mal positron fractions for E > 10 GeV are expected 
to be estimated with Ml and M2 models, respectively. 
We found that the MED and Ml models give similar 
positron fraction, so only the results with the MED and 
M2 models are shown in the following. In addition, we 
use b = l.Q x 10~ 16 x Eq cV GeV/sec. For the case where 
the primary positron is produced by the pair annihilation 
of dark matter, the source term is given by 
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Here, B is the so-called boost factor, (av) is the averaged 
pair annihilation cross section, and todm is the mass of 
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FIG. 1: Positron fraction for the Wino dark matter case with 
MED (top figure) and M2 (bottom figure) propagation mod- 
els. We take m w = 200 GeV (left) and 1 TeV (right), and 
the boost factor is taken to be B = 1, 3, and 10 (from bot- 
tom to top for each figure). The dotted-dashed line is the 
positron fraction calculated only with the background fluxes. 
The PAMELA data are also plotted. 



dark matter particle. In addition, phaio is the dark matter 
mass density for which we adopt the isothermal profile 



.11] p(r) = p (r c 2 ore + r|)/(r, 
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core 



+ *kpc) 



0.43 GeV/cm 3 is the local halo density, 



where ~ 
oro ~ 2.8 kpc 



is the core radius, r ~ 8.5 kpc is the distance between 
the galactic center and the solar system, and rk pc is the 
distance from the galactic center in units of kpc. If the 
positron is from the decay of dark matter, 
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where tdm is the lifetime of dark matter. In the above 
expressions, [dN e + / dE] ann and [dN e + /dE\& cc are the en- 
ergy distributions of the positron from single pair anni- 
hilation and decay processes, respectively, and are calcu- 
lated by using PYTHIA package [12| for each dark matter 
candidate. 

Now, we discuss the positron flux for several models 
of superparticle dark matter. We first consider the case 
where the LSP is Wino-like neutralino, which we denote 
W°. (Implications of the Wino dark matter scenario to 



the PAMELA results were first discussed in Q.) The 
thermal relic density of the Wino is much smaller than 
the present mass density of dark matter if < 1 TeV 
(with being the Wino mass). However, its mass den- 
sity can be consistent with the present dark matter den- 
sity if the Wino is non-thermally produced in the early 
universe [l3| (or if m^, ~ 2.5 — 3 TeV [14], which we 
do not consider in this letter). In addition, the Wino be- 
comes the LSP in some of the well-motivated scenarios of 
SUSY breaking, like the anomaly- mediation model [l5l ]. 
If W° is dark matter, it mainly annihilates into W + W~ 
pair when < 1 TeV. Then, the subsequent decays 
of W produce energetic positrons, which become the 
source of the cosmic-ray positron. 

We calculate the positron flux from such an annihi- 
lation process. In Fig. [1] we plot the positron fraction 
using the MED and M2 propagation models. (Results 
with the Ml model are similar to those with the MED 
model, and hence are not shown.) In the figure, the 
Wino mass is taken to be 200 GeV, and 1 TeV, for 
which (av) = 1.9 x 10" 24 , and 8.5 x 10~ 26 cm 3 /sec, 
respectively. In calculating the positron fraction, we 
adopt the following background fluxes [$ c -]bg = 
0.16^7(1 + 11E™ V + 3.2E%$) + 0.70£g c 7 v /(l + 
110^ v +600^-| v +580^v) GeV -1 cm" 2 sec" 1 str" 1 
for the electron, and [$ e +] B G = 4.5£&T V /(1 + 650.E^ V + 
1500£'qov) GeV -1 cm -2 sec -1 str -1 for the positron. As 
one can see, when the Wino-like LSP is dark matter, an 
enhancement of the positron flux is possible. However, 
irrespective of the boost factor, such a scenario fails to 
realize the increasing behavior of the positron fraction 
at 20 GeV < E < 100 GeV with the MED propaga- 
tion model. This is mainly because the positrons from 
hadrons produced by the hadronic decays of W ± signif- 
icantly contribute to the flux much below the thresh- 
old. On the contrary, with the M2 model which tends 
to suppress low-energy positrons relative to high energy 
ones, the shape of the positron fraction becomes consis- 
tent with the PAMELA data at high energy region (i.e., 
E > 20 GeV) for m w ~ 200 GeV and B ~ 3. How- 
ever, in such a case, the agreement between the theoret- 
ical positron fraction and the PAMELA data is poor at 
E < 20 GeV with the present choice of background. In 
addition, in the Wino dark matter case, the anti-proton 
flux tends to be too large [y] , which may exclude the pos- 
sibility of explaining the PAMELA positron excess in the 
Wino dark matter scenario. 

For a quantitative discussion, we calculate the x 2 vari- 
able using the PAMELA data [2j, neglecting effects of 
systematic errors in the theoretical calculation. Since 
the positron fraction in the low-energy region is sensitive 
to the background fluxes, we only use the data points 
with E > 15 GeV (5 data points) in the calculation 
of x 2 . (The detailed procedure to calculate x 2 is the 
same as those given in @.) With the M2 propagation 
model, the % 2 variable can become smaller than 11.0, 
which corresponds to the 95 % C.L. bound, in the pa- 
rameter region from = 130 GeV (with B ~ 1) to 
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FIG. 2: Positron fractions for the gravitino dark matter case. 
Gravitino is assumed to decay only into the first generation 
lepton (plus gauge or Higgs boson). We take m 3 / 2 = 250 GeV, 
500 GeV, and 1 TeV (from left to right), and the MED (top 
figure) and M2 (bottom figure) propagation models are used. 
For the case with the MED propagation model, we take r 3 / 2 = 
8.5 x I0 26 sec x (m 3 / 2 /100 GeV)" 1 , while for the case with 
the M2 model, we take t 3/2 = 1.8 x 10 2e sec, 9.0 x 10 25 sec, 
and 6.3 x 10 25 sec for m 3/2 = 250 GeV, 500 GeV, and 1 TeV, 
respectively. 

m w = 310 GeV (with B ~ 10). In addition, with MED 
and Ml models of propagation, the positron fraction ob- 
served by the PAMELA is hardly realized with Wino 
dark matter scenario with the conventional estimation of 
the background. 

Next, we discuss the case where the gravitino (denoted 
as VvJ i s the LSP and hence is dark matter. The pair 
annihilation cross section is negligibly small in such a 
case. However, with i?-parity violation (RPV), the grav- 
itino LSP becomes unstable and energetic positron can 
be produced by the decay. Even if the gravitino is un- 
stable, it can be dark matter if the RPV is weak enough 
so that the lifetime of the gravitino T3/9 is much longer 
than the present age of the universe [H, |T3|. In fact, 
such a scenario has several advantages. In the grav- 
itino LSP scenario with RPV, the thermal leptogenesis 
[HI becomes possible without conflicting the big-bang 
nucleosynthesis constraints. In addition, the fluxes of 



FIG. 3: Positron fractions for the gravitino dark matter case. 
Gravitino is assumed to decay only into the second generation 
lepton. We take m 3/2 = 250 GeV, 500 GeV, and 1 TeV (from 
left to right), and the MED (top figure) and M2 (bottom 
figure) propagation models are used. For the case with the 
MED (M2) propagation model, we take t 3 / 2 = 2.5 x 10 26 sec, 
1.6 x 10 26 sec, and 1.3 x 10 26 sec (1.6 x 10 26 sec, 1.0 x 10 26 sec, 
and 7.9 x 10 25 sec) for m 3/2 = 250 GeV, 500 GeV, and 1 TeV, 
respectively. 

the positron and 7-ray can be as large as the observed 
values, and the anomalies in those fluxes observed by 
the HEAT [|[ and the EGRET experiments, respec- 
tively, can be simultaneously explained in such a scenario 
if T3/2 ~ O(10 26 sec) [|[20j]. 

Here, let us consider the bi-linear RPV interactions. 
Using the bases where the mixing terms between the 
up-type Higgs and the lepton doublets are eliminated 
from the superpotential, the relevant RPV interactions 
are given by 

£rpv = B t .LiH u + m 2 iHd LiH^ + h.c, (4) 

where Li is left-handed slepton doublet in z-th gener- 
ation, while H u and Hd arc up- and down-type Higgs 
boson doublets, respectively. Then, the gravitino de- 
cays as W ~~ * 'fW^, ViZ, Vih, and ^7, where if 
and Vi are the charged lepton and the neutrino in i- 
th generation, respectively. Taking account of all the 
relevant Feynman diagrams, we calculate the branch- 
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ing ratios of these processes [9j. When the gravitino 
mass is larger than my/, the dominant decay 

mode is ipn ~ * In addition, in such a case, 

t 3/2 ~ 6xl0 25 secx(K l /10- 10 )- 2 (m 3/2 /l TeV)~ 3 , where 
Ki = (Bi sin (3 + m| cos 0) / to 2 . is the ratio of the vac- 
uum expectation value of the sneutrino field to that of 
the Higgs boson, with tan/3 = (H®)/(H%), and m„ i being 
the sneutrino mass. Thus, the lifetime of the gravitino is 
a free parameter and can be much longer than the present 
age of the universe if the RPV parameters Bi and to- „ 
are small enough. 

We calculate the positron flux from the decay of the 
gravitino dark matter. For simplicity, assuming a hier- 
archy among the RPV coupling constants, we consider 
the case where the gravitino decays selectively into the 
lepton in one of three generations (plus V4 7± , Z, or h). 

When the gravitino decays only into first genera- 
tion lepton, the dominant decay mode is ijj^ — ► e^W^ 
and a large amount of positrons with the energy of 
( TO 3/2 — rn w)/2 TO 3/2 are produced by the decay. Such 
monochromatic positron results in the significant increase 
of the positron fraction at 20 GeV < E < 100 GeV. 
Then, contrary to the Wino dark matter case, a dras- 
tic enhancement of the high energy positron fraction is 
possible with any of the propagation model. In such a 
case, the PAMELA anomaly is well explained with the 
MED model; we found that, with the M2 model, the 
high energy positron fraction is too much enhanced to 
be consistent with the PAMELA data. We found that 
the x 2 variable may become smaller than the 95 % C.L. 
bound (i.e., 11.0) with the MED model; if t 3 / 2 is prop- 
erly chosen, the positron fraction well agrees with the 
PAMELA data for to 3/2 > 100 GeV irrespective of the 
gravitino mass. (Simultaneously, the energetic 7-ray flux 
is also enhanced, which can be an explanation of the 
7-ray excess observed by the EGRET [H[.) With the 
MED model, we found that the the PAMELA data sug- 
gests t 3/2 ~ 8.5 x 10 26 sec x (to 3/2 /100 GeV)- 1 . In Fig. 
[21 the positron fraction is shown with this choice of the 
lifetime. We can see a good agreement between the the- 
oretical and observational positron fractions. Once the 
statistics will be increased, the PAMELA may see the 
end-point at E ~ (to 2 ^ 2 — to 2 v )/2to 3 / 2 . One can also 
see that the positron fraction at E ~ 10 GeV has a no- 
table dependence on the gravitino mass; with a better 
understanding of the background fluxes, it may be used 
to derive a more stringent constraint on the gravitino 
mass. 

When the gravitino decays into second or third gener- 
ation lepton, the positron from the decay of the primary 
lepton (i.e., /i + or t + ) is less energetic than the positron 
directly produced by the process tp^ — > e + W~ . Then, 
the energetic positron flux becomes suppressed. Even in 
such a case, it is still possible to have an enhanced high 
energy positron flux. For example, for the case where the 
gravitino decays only into the second generation lepton, 
we find the 95 % C.L. allowed region for any of the prop- 
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FIG. 4: Positron fractions for the case where dark matter de- 
cays only into e + e~ pair. We take ttidm = 250 GeV, 500 GeV, 
and 1 TeV (from left to right), and r DM = 2.2 x 10 27 sec x 
(m DM /100 GeV)" 1 . The MED propagation model is used. 

agation model. For the M2 model, which gives a better 
fit to the data than MED and Ml, the best-fit value of 
t 3/2 is 1.6 x 10 26 sec, 1.0 x 10 26 sec, and 7.9 x 10 25 sec 
for to 3/2 = 250 GeV, 500 GeV, and 1 TeV, respectively. 
In Fig. [31 we show the positron fraction for these cases. 
With the MED propagation model, the agreement be- 
tween the theoretical prediction and the PAMELA data 
becomes slightly worse. 

In the scenario with ip^ — * Z ± W =F , cosmic-ray anti- 
proton is also produced. However, we have checked that 
the resultant anti-proton flux becomes comparable or 
smaller than the observed values, taking account of the 
uncertainties in the propagation model of anti-proton. In 
particular, the best-fit lifetime to explain the PAMELA 
anomaly is about 5 times longer than that used in the 
calculation of 20], and we obtain the anti-proton flux 
smaller than that presented in [2(j. The detailed analy- 
sis is given in [2l| . 

So far, we have concentrated on two important candi- 
dates for SUSY dark matter. Even in other cases, how- 
ever, the positron flux may be also enhanced. For ex- 
ample, with the RPV interaction given in Eq. (J4|, the 
PAMELA anomaly may be explained for the Bino dark 
matter case. In such a scenario, the Bino dominantly 
decays as B — > e ± W T and the monochromatic positron 
produced by the decay becomes the origin of the ener- 
getic cosmic-ray positron. 

Our study indicates that the PAMELA anomaly may 
be well explained if (almost) monochromatic positrons 
are emitted from the decay or the pair annihilation of 
dark matter. Thus, we finally consider the case where 
the dark matter is unstable and dominantly decays into 
e + e~ pair. This is the case if the left- or right-handed 
sneutrino is the LSP and hence is dark matter, and also 
if the LiL\E\ type RPV superpotential exists. (The 
left- and right-handed sneutrinos are viable candidates 
for dark matter; see [13, HH-) With the \ 2 analysis, 
we found that the positron fraction becomes consistent 
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with the PAMELA data when r DM — 2.2 x 10 27 sec x 
(?ti dm /100 GeV)- 1 for m DM > 100 GeV in the MED 
model. In Fig. 01 we show the positron fraction for 
m D M = 250 GeV, 500 GeV, and 1 TeV, using the 
above tdm- As one can see, in such a case, an excel- 
lent agreement between the theoretical prediction and 
the PAMELA data is obtained. 
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